Universality Class of Thermally Diluted Ising Systems at Criticality 
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The universality class of thermally diluted Ising systems, in which the realization of the disposition 
of magnetic atoms and vacancies is taken from the local distribution of spins in the pure original Ising 
model at criticality, is investigated by finite size scaling techniques using the Monte Carlo method. 
We find that the critical temperature, the critical exponents and therefore the universality class 
of these thermally diluted Ising systems depart markedly from the ones of short range correlated 
disordered systems. Our results agree fairly well with theoretical predictions previously made by 
Weinrib and Halperin for systems with long range correlated disorder. 



I. Introduction 

During last decades the systems with quenched ran- 
domness have been intensively studied [Q. The Harris 
criterion |Q predicts that weak dilution does not change 
the character of the critical behavior near second order 
phase transitions for systems of dimension d with spe- 
cific heat exponent lower than zero in the pure case (the 
so called P systems), apure < =4> Vpure > 2/d, being 
ly the correlation length critical exponent. This crite- 
rion has been supported by renormalization group (RG) 
and scahng analyses d). The effect of strong di- 
lution was studied by Chayes et al. 0. For apure > 
(the so called R systems), the Ising 3D case for exam- 
ple, the system fixed point fiows from a pure (undiluted) 
fixed point towards a new stable fixed point §|,|,|,@,| 
at which arandom < 0. Recently Ballesteros et al. have 
used the Monte Carlo approach to study the diluted Ising 
systems in two three and four dimensions [ pd| . 
The existence of a new universality class for the random 
diluted Ising system (RDIS), different from that of the 
pure Ising model and independent of the average density 
of occupied spin states (p), is proved, using an infinite 
volume extrapolation technique pO| based upon the lead- 
ing correction to scaling. The critical exponents obtained 
this way could be compared with the experimental crit- 
ical exponents for a random disposition of vacancies in 
diluted magnetic systems Q. 

In all cases previously mentioned frozen disorder was 
always produced in a random way, that is, vacancies were 
distributed throughout the lattice randomly. Real sys- 
tems, however, can be realized with other kinds of dis- 
order, where the vacancy locations are correlated. In 



particular, long range correlation (LRC) has been found 
in X-ray and Neutron Critical Scattering experiments in 
systems undergoing magnetic and structural phase tran- 
sitions • This effect has been modeled by assuming 
a spatial distribution of critical temperatures obeying a 
power law g{x) ^ for large separations x p5[ |. In 
general these systems behave in a way very different from 
RDIS, since systems with randomly distributed impuri- 
ties may be considered as the limit case of short range 
correlated (SRC) distributions of the vacancies. The ba- 
sic approach to the critical phenomena of LRC systems 
was established by Weinrib and Halperin ||l^ almost two 
decades ago. They found that the Harris criterion can be 
extended for these cases, showing that for a < d, the dis- 
order is irrelevant if aVpure — 2 > 0, and that in the case 
of relevant disorder a new universality class (and a new 
fixed point) with correlation length exponent v — 2/a, 
and a specific heat exponent a = 2{a — d)/a appears. In 
contrast, if a > d, the usual Harris criterion for SRC sys- 
tems is recovered. LRC disorder has been studied also 
by the Monte Carlo approach . In this case a correla- 
tion function g{x) = a;~° with a = 2 (defects consisting 
in randomly oriented lines of magnetic vacancies inside a 
three dimensional Ising system) confirmed the theoretical 
predictions of Weinrib and Halperin. 

In the present paper we will study Ising three dimen- 
sional systems where the long range correlated dilution 
has been introduced as a thermal order-disorder distri- 
bution of vacancies in equilibrium, governed by a char- 
acteristic ordering temperature (6*), in a way similar to 
the thermal disordering in a binary alloy of magnetic 
(spins) and non-magnetic atoms (vacancies). [See pre- 
vious work fl^ and also the application of this kind of 
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disorder in percolation problems [Q]. Making more ex- 
plicit the analogy with order-disorder in alloys we may 
distinguish clearly between: 

i) thermally diluted Ising system realizations 
(TDIS), in which the quenched randomness is produced 
by considering a ferromagnetic Ising system ai {6 = 
T^^): after thermalization, the spins of the dominant 
type (concentration c > 0.5) are taken as the locations of 
the magnetic atoms, and the rest are taken as the mag- 
netic vacancies. The structure of the realization is fixed 
thereafter for all temperatures at which the magnetic in- 
teractions are subsequently investigated. 

it) random diluted Ising system realizations (RDIS) 
(equivalent to >> Tj?^), also fixed thereafter for all 
temperatures at which the spin interactions are inves- 
tigated. Of course, we chose for comparison a vacancy 
probability p = 0.5, resulting in c w 0.5. 

iii) as a third possibility, not considered here, one could 
investigate what might be called an antiferromagnet- 
ically diluted Ising system (with 9 << , which 
would lead to a disposition of non-magnetic atoms (va- 
cancies) strictly alternating with spins, with c = 0.5. 

So, if this ordering temperature determining the par- 
ticular realization is high enough, the equilibrium ther- 
mal disorder will be very similar to the random (short 
range correlated) disorder of the usual previous inves- 
tigations. On the other hand, if 6 happens to coin- 
cide with the characteristic magnetic critical tempera- 
ture {Tj}^ — 4.511617) of the undiluted system, we will 
have vacancies in randomly located points, but with a 
long range correlated distribution. [Note that the situa- 
tion differs markedly from that of previously studied LRC 
systems, in which lines or planes of vacancies were con- 
sidered] . The correlation of our TDIS is given by a value 
a = 2 — rjpure , where rjpure is the correlation function ex- 
ponent for the pure system. Since d = 3 and rjpure = 0.03 
for the three dimensional Ising system, we have a long 
range correlated disorder with a — 1.97 < 3 = d. So we 
are in the case where LRC disorder is relevant and we 
should detect a change of universality class with respect 
to the SRC case (following Weinrib and Halperin we ex- 
pect for the thermally diluted Ising system v ^ 1 and 
a « —1). Details about the construction of these ther- 
mally diluted Ising systems (TDIS) can be found in Ref. 

In the present work we study the critical behav- 
ior and the university class of three dimensional TDIS at 
criticality using the Monte Carlo approach. We will com- 
pare our results with the critical behavior of the RDIS. 

The structure of the paper is as follows: In Section 
II we study the dependence of the critical temperature 
(and of the self-averaging at criticality) with the size of 
the system for both TDIS and RDIS. Once the critical 
point is determined, we investigate whether or not TDIS 
and RDIS belong to the same universality class. In order 
to proceed, we study the critical behavior of both kinds of 
systems by applying finite size scaling techniques (Section 



III) and by using the effective-exponent approach (Sec- 
tion IV). A summary of the main results and concluding 
remarks are given in Section V. 

II. Transition temperature and self-averaging of 
thermally diluted systems 

For a hypercubic sample of linear dimension L and 
number of sites N = L"^, any observable singular property 
X has different values for the different random realiza- 
tions of the disorder, corresponding to the same dilution 
probability p (grand-canonical constraint). This means 
that X behaves as a stochastic variable with average X 
(in the following, the bar indicates average over subse- 
quent realizations of the dilution and the brackets indi- 
cate MC average). The variance would then be (AX)^, 
and the normalized square width, correspondingly: 



Rx = {AXf/X^ (1) 



A system is said to exhibit self-averaging (SA) if Rx 
as L — > oo. If the system is away from criticality, L » ^ 
(being ^ the correlation length). The central limit the- 
orem indicates that strong SA must be expected in this 
case. However, the behavior of a ferromagnet at critical- 
ity (with ^ >> L) is not so obvious. This point has been 
studied recently for short range correlated quenched dis- 
order. Aharony and Harris (AH), using a renormalization 
group analysis in d = 4 — e dimensions, proved the expec- 
tation of a rigorous absence of self-averaging in critically 
random ferromagnets |2^]. More recently, Monte Carlo 
simulations were used to investigate the self-averaging in 
critically disordered magnetic systems |l0| , ^ , ^ . The 
absence of self-averaging was confirmed. The normalized 
square width Rx is an universal quantity affected just 
by correction to scaling terms. LRC diluted systems are 
expected to have different critical exponents and differ- 
ent normalized square widths with respect to those of the 
usual randomly disordered systems studied previously. 

We perform Monte Carlo calculations of the magne- 
tization and the susceptibility (x = {(M^) - {Mf)/T) 
per spin at different temperatures for different realiza- 
tions of TDIS, and for randomly diluted systems with 
p ~ 0.5 (restricting to c > 0.5) using in both cases 
the Wolff ||2^ single cluster algorithm ||2^ with peri- 
odic boundary conditions, on lattices of different sizes 
L = 10,20,40,60,80,100. Resuhs for susceptibility vs. 
temperature are shown in Fig. 1. 
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FIG. 2. Semi log representation of critical temperature Tc 
vs. lenght L for random dilution (p = 0.5) (black) and for 
thermal dilution (white). The continuous line is the fit ob- 
tained using 1/ = 0.6837 (short range correlated random expo- 
nent) and dashed line the fit obtained using u = 1.015 (long 
range correlated exponent, a=1.97) in either case. 

In both cases a fit to Eq. 2 has been performed for 
both values v — Vrandom — 0.683 (continuous line) and 

— ^thermal — 1.015 (dashed line). Note how the ther- 
mal data fit better Eq. 2 for Tc{L) using vthermah indi- 
cating a possible change in universaUty class with respect 
to the random case. [If we fit the data leaving all param- 
eters free, we find lyrandom ~ -0.7 and vthermai ~ 1-2, 
which are very near the expected results]. The extrapo- 
lated values of critical temperatures for infinite systems 
obtained this way are TJ''"'''""(oo) = 1.845 ±0.003 (close 
to the values previously obtained by Ballesteros et al. 
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FIG. 4. Scattered plot of magnetization M at the criti- 
cal temperature vs. concentration c for the realizations con- 
sidered of the random (p = 0.5) and the thermal dilutions. 
(L = 20,40,60,80). 

From Fig. 4 we can extract averaged values for the 
magnetization and the inverse susceptibility, M and 
(with X = at the critical point) for the TDIS. Both 

averaged values are expected to fit the following scaling 
laws at criticality: 

mV) cx L-^'" (3) 

FHiy^L-^/^ (4) 

Considering {l/v)thermai — 0-12 — 0.985, we can ob- 
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(c(L) -0.5) oc L-('3/'^)3o 



(6) 



□ 



□ 



where {P/ 1^)30 — 0.52 gives the values corresponding to 
the pure three dimensional Ising case, because in criti- 
cally thermally diluted Ising systems, vacancies are dis- 
tributed with the same long range correlation spin 
distribution function as in the pure case The fit to 
the average value of the concentrations shown in Fig. 4 
give, for the thermal value {P/ 1^)30 ^ 0.55 ±0.08. 

This implies also a clear difference between RDIS and 
TDIS, since Ising systems with vacancies randomly dis- 
tributed are not expected to follow an scaling behav- 
ior with {/3/iy)3D- [Fitting to an scaling law the results 
for RDIS gives an exponent around 1.4, which implies a 
much faster convergence to c = 0.5]. 



IV. Effective exponents and possible crossover 
for finite systems 

The difference between the universality class of RDIS 
and TDIS can be detected also by means of the effective 
critical exponents. In the case of the magnetization the 
effective critical exponent is defined by: 



^ d log(A/)/a log(0 - - T) 



(7) 



X 



FIG. 5. Average value of the inverse of the susceptibility 
~^{L) vs. the lenght of the system L for the random dilution 



(p = 0.5) case (black) and the thermal dilution case (white). 
Continuous line and dashed line indicate fits to the random 
and the thermal cases, respectively. 

The critical exponents obtained are in agreement with 
those of Ballesteros et al. within our error bars. 

The average concentration for the thermal system is 
also expected to show a scaling law behavior given by: 



with the critical temperature of the particular real- 
ization (i) characterized by a maximum of the suscep- 
tibility (x = (m2) - (Mf ). For L ^ 00 and t 0, 
(3eff — P- Finite size effects force the effective critical 
exponent to drop to zero before the critical value is at- 
tained. However, since Prandom — 0.3546 (Ref. pO|] ), 
and Pthermai IS cxpcctcd to be around 0.5 (Ref. |16| ) 
, effective critical exponents (for the thermal case) may 
rise to values greater than 0.3546 and lower than 0.5, 
before finite size effects appear, indicating the difference 
of universality class between both kinds of systems. In 
the RDIS the effective critical exponent is expected to 
be always lower than 0.3546. Monte Carlo simulations 
of magnetization vs. temperature have been performed 
for randomly [p ~ 0.5) and thermally diluted systems, 
with L = 80. In Fig. 6 we show the results for P^ff 
vs. log(t) for two samples of the TDIS and the RDIS 
type respectively. The random effective critical expo- 
nent is always under 0.35 and it seems to tend towards 
this value for large enough L, as expected, but for ther- 
mal systems the behavior is completely different: In the 
figure the value of the critical thermal effective exponent 
is between 0.35 and 0.5. An analogous investigation has 
been done for different values of L and different realiza- 
tions. The same effect has been found for lower L val- 
ues. However, we may note that in TDIS the effective 
critical exponent arrives at the maximum in a very dif- 
ferent way depending on the realization. The reason is 
twofold: (1) the different disposition of the vacancies in 
each particular realization, and (2) the large differences 



5 



the size of the fixed vacancies and spin clusters, the size 
of the system itself and the thermal spin fluctuations 
might all play a role. The possibility that for L — > oo 
and c = 0.5 the apparent crossover could disappear all 
together should not be excluded, entailing no crossover 
from critical thermal to random universality class. 

V. Concluding remarks 

To summarize, a new way to produce diluted Ising sys- 
tems with a long range correlated distribution of vacan- 
cies has been analyzed by means of Monte Carlo cal- 
culations and finite size scaling techniques. We find an 
universality class different from that reported for diluted 
Ising systems with short range correlated disorder. Our 
systems may be included in the universality class pre- 
dicted by Weinrib and Halperin for a « 2. This kind 
of thermal disorder had been already applied in percola- 
tion problems, but as far as we know it had never been 
applied to magnetic systems, in which long range corre- 
lated disorder has been previously produced mostly by 
random distribution of lines or planes of vacancies. The 
present dilution procedure may be applicable to systems 
where the long range correlated disorder is not due to 
dislocations, preferential lines or planes of vacancies, but 
to systems where the vacancies (points) are critically dis- 
tributed in clusters as in the case of order-disorder sys- 
tems. 
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erous access to his computing facilities. We acknowl- 
edge financial support from DGCyT through grant PB96- 
0037 and from the Basque Regional Government (J.I.). 



FIG. 6. Effective critical magnetization exponent /3e// vs. 
logit), with t = Tl — T, for the thermal dilution case (wliito) 
and the random dilution (p = 0.5) case (black). The realiza- 
tions considered correspond to L = 80. 

In Fig. 6 the effective critical exponent for the thermal 
case seems to produce a crossover towards /3 = 0.35 as the 
temperature approaches the critical point, before finite 
size effects finally appear. In principle this might be influ- 
enced by the indetermination in the measurement of the 
particular critical temperature of the realization. How- 
ever, this crossover may point out that different length 
scales might be important at different distance from the 
critical point. In principle characteristic lengths such as 
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